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ABSTRACT 

Recent radio pulsar observations have shown that a number of pulsars display interest- 
ing long term periodicities in their spin-down rates. At least some of these pulsars also 
undergo sharp changes in pulse profile. This has been convincingly attributed to the 
stars abruptly switching between two different magnetospheric states. The sharpness 
of these transitions has been taken as evidence against free precession as the mecha- 
nism behind the long term variations. We argue that such a conclusion is premature. 
By performing a simple best-fit analysis to the data, we show that the relationship 
between the observed spin and modulation periods is of approximately the correct 
form to be accounted for by the free precession of a population of neutron stars with 
strained crusts, the level of strain being similar in all of the stars, and consistent 
with the star retaining a memory of a former faster rotation rate. We also provide 
an argument as to why abrupt magnetospheric changes can occur in precessing stars, 
and how such changes would serve to magnify the effect of precession in the timing 
data, making the observation of the precession more likely in those stars where such 
switching occurs. We describe how future observations could further test the precession 
hypothesis advanced here. 
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1 INTRODUCTION 

Radio pulsars are remarkably stable rotators. However, as 
was apparent from soon after the discovery of the first pul- 
sar, in addition to a smooth secular spin-down, they exhibit 
other interesting timing features, including sudden spin-ups 
known as glitches, periodic variations in spin-down rate pos- 
sibly attributable to free pression, and low-level random tim- 
ing fluctuations known as timing noise. 

The study of all these forms of spin-down variation are 
of great interest as probes of neutron star structure. The 
glitch and precession phenomena are believed to depend sen- 
sitively on the coupling bet ween the s o lid an d superfluid 
phases of the star (see e.g. iLink et al.l (|1999r 0. while the 
origin of the timing noise is less well understood, and may 
be caused by process interna l to the star or external (i .e. 
in the pulsar magnetosphere IjCordes fc Greensteinlll98ll )). 
The reported existence of sustained long-period free preces- 
sion is of particular interest, as such motions are inconsistent 
with the pinned superfluid model that has become the stan- 
dard explanation of large glitches l| Jones fc A ndcrsson 2001; 
ILink fc Epsteinll200ll : lLinkll2003l . 120061) . Additionally, timing 
noise is expected to limit the performance of pulsar timing 
arrays, whose aim is to detect low fre quency (nano-Hertz) 
gravitational waves (|Jenet et al.ll2005T ). As pointed out re- 



cently, if timing noise could be modelled and removed from 
the timing solution, the sensitivity o f pulsar timing arrays 
could be improved (ILvne et alj 120101 ). There are therefore 
both theoretical and practical motivations for understand- 
ing better the origin of timing noise. 

As summarised below, recent radio observations have 
shown that some pulsars with harmonic features in their 
timing residuals undergo sudden sharp changes in pulse pro- 
file which, in one case at least, correlates with sharp change s 
in spin-down torque IjKramer et alj|200a : ILvne et alj|20ich . 
As was convincingly argued, this is good evidence for the 
star switching back and forth between two different magne- 
tospheric states. The sharpness of this switching has been 
taken as evidence that free precession is n ot a viable mecha - 
nism to explain the harmonic timing data IjLvne et alj201Ch . 
thereby removing the theoretical difficulty discussed above. 

However, as we argue below, we believe it is premature 
to abandon the precession hypothesis as an explanation for 
the harmonic timing data. We give two main arguments to 
support this. Firstly, the ratios P/P mo d of the spin periods 
P (of order a second) to the long term periodicities P mo d 
(of order a few years) are of the right magnitude to be ac- 
counted for by free precession of a star whose crustal stress 
is due to it retaining a memory of a former faster rotation 
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rate; few (if any) other known processes act on such long 
timescales. We find that the implied breaking strains must 
be at least as large as 5 x 10 _ , with minimum birth spin 
frequencies of 50 Hz if only the crust participates in the pre- 
cession; these numbers increase to 5 x 10 -2 and 500 Hz if 
the whole star precesses as one. That free precession of a 
strained star, with these sorts of minimum strains and birth 
frequencies, is a viable explanation for the timing behaviour 
of PSR B1828-11 has already been demonstrated in detail 
bv lCutler et al.l (|2003r >: what we add here is the demonstra- 
tion that the scaling of P mo d with respect to P extends over 
a whole sample of potentially precessing pulsars, consistent 
with a similar level of strain in each. 

Secondly, we argue the sharp magnetospheric switching 
and precession are not mutually exclusive. Our argument 
is qualitative, but, we believe, persuasive. In brief, suppose 
some pulsar is delicately balanced between two magneto- 
spheric states. This would be the case if the energy available 
to accelerate particles is close to some threshold for some 
avalanche-like process, say pair production. Then, if the star 
is precessing, at some precessional phases the critical energy 
for the process is exceeded, favouring one magnetospheric 
state, while at other precessional phases the critical energy 
is not exceeded, favouring the other magnetospheric state. 
That some pulsars should be so delicately balanced is plau- 
sible, given the random nulling and mode changing observed 
even in non-precessing pulsars. The precession simply adds 
a bias to the probability of the magnetosphere being in one 
state or the other. The corresponding changes in spin-down 
torque would then serve to amplify the effect of precession on 
the timing data: it would then be no surprise that switching 
behaviour and precession have been observed to go hand-in- 
hand. 

It is this precession hypothesis that we explore in this 
paper, whose structure is as follows. In section [5] we give a 
summary of the relevant observations. In section[3]we set out 
our model. In section [4] we describe some consistency tests 
that we have applied to our model. In section [S] we discuss 
our results, and comment on some alternative explanations. 
Finally, in section [6] we summarise our findings and give 
some concluding remarks. Where relevant, we have assumed 
neutron star masses of M = 1.4Mq, radii R — 10 6 cm, and 
moments of inertia / = 10 45 g cm 2 . We have made use of 
the AT NF pulsar database, as described in lManchester et al.l 
(20050 



radio pulsations, or off, contributing to neither spin-down 
nor radio pulsation (|Kramer et al.ll200q ). 



More recently, iHobbs et al.l (J2010j) have presented the 
results of a study of timing noise in 366 pulsars, using data 
sets gathered at Jodrell Bank. The long duration of the data 
sets (~ several decades) enabled them to identify pulsars 
whose timing residuals displayed interesting long-period (~ 
several year) har monic structure. T his theme was taken up 
in more detail bv lLvne et al.l (J2010J). The following interest- 
ing features emerged: (i) A subset of 17 pulsars displayed sig- 
ni ficant harmonic va riation in their spin-down rate (see Fig 2 
of lLvne et al.l (120101 0. (ii) In at least six of these pulsars, the 
variations in the spin-down rate were correlated w ith varia- 
tions in pulse profile (see Fig 4 of lLvne et al.l (J2010J)). On the 
basis of these two observations, one might be tempted to ad- 
vance free precession as the mechanism responsible for the 
variations. Howeve r, motivated by th e peculiar properties 
of PSR B1931+24, iLvne et all (|201Ch analysed the pulsars 
further and found (iii) For at least two pulsars, the varia- 
tions in pulse profile were smooth only when averaged over 
a sufficiently long timescale (~ 100 days) . When studied on 
shorter timescales (~ 10 minutes), the pulse pr ofile showed 
signs of being in one of two states (see Fig 2 of lLvne et al.l 
1)20101) '). with the switch from one state t o another occurrin g 
on a very short timescale; see Fig 5 of ILvne et al.l IJ2010T) . 
The longer-term smooth harmonic variation was then to be 
understood as being caused by the fraction of time the star 
spends in one state or the other smoothly varying on the 
long modulation timescale. 

This last point, viz the sharp switching between pulse 
profiles, is clearly reminiscent of the behaviour of PSR 
B1931+24. This in turn suggested that not only the pulse 
profile but also the spin-down torque itself in the set of 
17 pulsars was undergoing sudden sharp changes, caused 
by a similar sort of magnetospheric switching process as 
in B1931+24, although it was not possible to confirm this 
as torque variations can be detected only via variations in 
spin-down rate, and the spin-down could only be measured 
by averaging over timescales much longer than the switch- 
ing timescale. In any case, the precessional hypothesis, with 
its smoothly varying modulation in pulse profile and spin 
parameters (see below) was rejected in favour of a magneto- 
spheric switching hypothesis, which more naturally seems to 
accom modate sharp chan ges in spin-down torque and pulse 
profile (|Lvne et alj|2010h . 



2 SUMMARY OF RELEVANT OBSERVATIONS 



iKramer et all (120061 ) have reported on PSR B1931+24, 
which has the unique property of displaying large and cor- 
related changes in spin-down rate and radio pulse strength: 
the pulsations are either 'on' or 'off', and the star seems to 
spin-down approximately 50% faster when on than when off. 
The switching between states occurs sharply-on a timescale 
of less than 10 seconds. This has a natural interpretation 
as being caused by some magnetospheric process, which is 
either on, contributing to the spin-down and providing the 



http : //www . atnf . csiro . au/resear ch/pulsar/psrcat 



However, it remains the case that long term periodic- 
iti es are seen i n the spin-down of the 17 pulsars discussed 
bv lLvne et al.l (|2010r ). and that some have such clear timing 
variations that one can easily pick them out by eye, e.g. see 
the plots of spin-down rate for PSRs B 1540-06, B1642-03 , 
B1826-17 and B1828-11 in Figure 2 of ILvne et ail i|2010h . 
There is no obvious mechanism to provide the clock for this 
within the magnetospheric switching hypothesis as formu- 
lated above. Given this deficiency, we will now reexamine 
the pre cession hypo t hesis as an explanation for the timing 
data of ILvne et al.l 1)20101 ). paying particular attention to 
whether or not free precession is likely to be consistent with 
the sharp switching in pulse amplitude seen in some of the 
pulsars. 
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3 BASIC MODEL 

Let us suppose that free precession is indeed the source of 
the harmonic structure seen in the timing data. That free 
precession should produce variations in pulse arrival time is 
well known l]Rudermanl . 1970b; Goldreich 1970). Following 
iGoldreichl |l970j), if a neutron star is (for the sake of sim- 
plicity) idealised as a biaxial body with principal moments 
of inertia [7i,7i,73], its angular velocity vector Oi evolves 
according to 

(1) 
(2) 
(3) 

where Ti is the spin-down torque. The motion of the body 
can then be understood using a method of successive approx- 
imations. To leading order the motion is the well known free 
precession, found by setting T = in the above. Assum- 
ing that the body is close to spherical ((73 — 7i)/7i <C 1) 
and that the wobble angle 8 (the constant angle between 
the body's 3-axis and its fixed angular momentum vector) 
is small (fl<l) one obtains 

fii = oe cos ip, (4) 

2 = 00 sin V, (5) 

n 3 = n, (6) 

where O = (f2j£2j) ' 2 is the magnitude of fij, ip is a phase 
that tracks the slow precession motion 



ip — ftfpt, 
Clip the slow precession frequency 

fifp = £31 Cl, 



(7) 



(8) 



£3i the fractional asymmetry in the moment of inertia tensor 
7 3 - 7i 



631 



<L 



(9) 



and 7 denotes the stellar moment of inertia, neglecting the 
small difference between different axes. 

This torque-free precession produces variations in the 
timing and pulse profiles. The apparent spin frequency is 
modulated by a fractional amount £31 9 once per precession 
cycle Pip = 2ir/Clfp, while the pulse shape will be modulated 
by the pulse emission cone wandering over a latitudinal an- 
gular interval ±8 about its aver age location, relative to a 
fixed observer (|Rudermanll 19700 ). 

This torque-free motion can then be inserted into the 
Ti 7^ equations of motion given above to compute the 
motion to the next order of accuracy. The variation in spin 
frequency is then 
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l |GoldreicrJll970h . As was emphasised bv lCordej l|l993h . the 
torque, in addition to producing the well-known secular spin- 
down, introduces additional timing variations: the free pre- 
cession modulates the spin rate and the latitudinal position 
of the magnetic axis. Assuming that the spin-down torque is 
sensitive to such changes, this produces a variability in the 
torque, i.e. a time variation in the torque locked in phase 
with the precession, which in turn modulates the spin-down 
rate. The precise form of the spin-down variation depends 



upon the functional form of T. Schematically, we can ne- 
glect the small last term of equation (|10|) to give the (rather 
obvious) relation 
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where 8T[<j>(t)] denotes that part of the torque that is 
modulated by the precessional motion. We have written 
ST = 8T[(f)(t)\ to make it clear that the modulation is a func- 
tion of precessional phase, which is itself a linear function of 
time. Integrating with respect to time gives the correspond- 
ing variation in spin frequency 



SO(t) = 



5T[<t>(t)] 



dt 



and integrating again gives the phase residuals: 
c5$(i) = f SO(t) di. 



(12) 



(13) 



If the variation with torque is harmonic, i.e. if ST is 
an oscilliatory function of time with frequency 0{ p , these 
integrals with respect to time simply introduce factors of 
the free precession period fif p 

ST - <)T •• ' )T (14) 
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These equations can be rewritten in a more illuminating 
form. Defining a spin-down timescale r s( j = fi/(2fi), and 
writing in terms of the fractional variation in torque ST jT 
and free precession period Pf p — 2n/Oi p , we obtain 
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The torque itself is believed to be composed of two 
parts: a part from direct emission of electromagnetic waves 
from the time-varying magnetic multipole moments of the 
star, and a part caused by electromagnetic emission from 
charged particles being accelerated in the magnetosphcrc 
(Goldreich & Julian 1969; Ruderman & Sutherland 1975). 

The direct electromagnetic emission from the star's 
magnetic multipoles is probably dominated by the magnetic 
dipole, explicit forms of which have b een worked out previ- 
ously (|Goldreichlll970l ; iMelatos 20001). In any case, for each 
contributing multipole, the associated spin-down torque will 
be some function of the time derivatives of a 'magnetic axis' 
rrii, fixed in the star, and the time derivatives themselves 
are generated by the rotation Ot: 

d n+1 rrn d n m k 

1 ^r = ^n i -^ r . (16) 

The precessional motion will smoothly modulate these time 
derivatives (see equations ©-(O), so that this piece of the 
spin-down torque will be a smooth function of the preces- 
sional phase and therefore a smooth function of time. We 
can then make a simple estimate ST/T ~ 9, leading to 
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It was equations of this form, adjusted to reflect the exact 
functional form of the magnetic dipole sp in-down torque, 
that were used bv ljones fc Anderssonl (J200l|) to extract wob- 
ble angles 9 from the timing variations of precession candi- 
dates. 
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Now let's consider the magnetospheric contribution to 
the spin-down torque. Unfortunately, the magnetospheric 
torque is based on intrinsically more complicated physics 
than the direct electromagnetic emission, so explicit formu- 
lae for the torque as a function of stellar magnetisation and 
spin rate are not available. Nevertheless, it is straightfor- 
ward to argue that the magnetospheric torque should have 
a part locked in phase with the precession. Ultimately, the 
magnetospheric torque requires the extraction of charged 
particles from the stellar surface by the strong electric fields 
Ei generated by the rotation of th e magnetic field Bj embed- 
ded in the highly conducting st ar l|Goldreich fc Julianl! 19691 : 
iRuderman fc Sutherland! 1 197a ). Assuming perfect conduc- 
tivity at the stellar surface, the induced electric field is 



E 



v x B, 



(18) 



where v = CI x r is the velocity of the stellar surface, i.e 
1 



E 



-(SI x r) x B. 



(19) 



Given that, to an excellent approximation, Bi will be fixed 
with respect to the star, the dependence of fli on preces- 
sional phase guarantees the precession-sensitive nature of 
this accelerating electric field. This modulation in the accel- 
erating field will inevitably produce modulation in the ener- 
gies to which particles are accelerated in the magnetosphere, 
and therefore also to modulation in the spin-down torque, al- 
though the link between these depends upon the exact and 
poorly understood magnetospheric mechanism at work. It 
follows that even in the case of magnetospheric torques, the 
spin-down rate will be modulated by the precession motion. 

As described in section [2] there is evidence for sharp 
switching behavior in the pulsar profile, which has been 
argued to correspond to sharp changes in the spin-down 
torque. How can this be reconciled with free precession, 
where the spin variation of equations ((3HS1) above are smooth 
functions of time? In the case of spin-down by direct elec- 
tromagnetic emission, there is clearly no way of producing 
sharp changes in torque: any magnetic multipole tied to the 
star's precessional motion will have a smooth time variation 
(as viewed from the inertial frame) and so will produce a 
smooth variation in torque. In the case of magnetospheric 
torque, the situation is less clear. One might be tempted to 
argue that the observed switching represents the precession 
pushing the available particle energy at some point in the 
magnetosphere over some critical threshold for some emis- 
sion process to occur. For those precessional phases where 
the particle energy exceeds the threshold, the process oc- 
curs, while for those precessional phases where the energy is 
too low, the process cannot occur. 

However, this in itself will not lead to sharp switch- 
ing, as presumably the area of stellar surface over which the 
threshold is exceeded is itself a smoothly varying function 
of precession phase, going from zero up to some maximum 
and smoothly back to zero once per precession cycle. What 
is needed for the sharp switching is that, once the energy 
threshold is exceeded at a single point, the process rises to 
a significant level, i.e. there is an 'avalanche' effect. Pair 
production is a clear candidate for providing this avalanche 
behaviour. The energy threshold would be that for pair pro- 
duction itself, related to the non-zero rest mass of the elec- 
tron/position pair, although there are other energy thresh- 



olds that might be relevant, such as the work function for 
particle extraction from the stellar surface. 

That the magnetosphere might be so delicately bal- 
anced is made plausible by t he well known phenomenon o f 
mode changing and nulling (ILvne fc Graham-Smith! 120061 ). 
These are well established phenomenon, showing that even 
in non-precessing stars the magnetosphere can be delicately 
balanced between a few quasi-equilibrium states, with ran- 
dom changes between the states occurring. The addition of 
precession will then presumably weight the magnetosphere 
toward one state or another, depending upon the preces- 
sional phase. That the switching in precessing stars is only 
quasi-randomly related to the precession phase is then a 
natural consequence of the random nature of mode switches 
in non-precessing stars. Our proposal is that the precession 
then causes the star to move between phases where it is more 
likely (as compared to average) to be in one magnetospheric 
state and phases where it is less likely to be in the state. 

In our precession model this corresponds to a compo- 
nent ST of the torque being either on or off, with a probabil- 
ity tied to the precessional phase. In this case we no longer 
have ST ~ 9T; the size of the torque variation is determined 
not by the precession angle 9 but by the magnitude of the 
process being switched on/off in the magnetosphere. The 
wobble angle only needs to be large enough to carry the star 
from below the threshold to over the threshold (and back) 
once per precession cycle. It is still the case that equations 
(|ll[l - (|13p apply, but now ST is some function only loosely 
tied to the precessional phase tp. That such a fluctuating 
to rque can account for real timing data was demonstrated 
bv lLvne et all |20ld ) for PSR B1828-11; see their Figure SI. 
Given that this is an additional source of harmonic varia- 
tion, over and above the smooth variation, this will tend to 
make precession more noticeable in pulsars with switching 
magnetospheric states. 

To sum up, our model is as follows. The basic mech- 
anism at work in creating the periodicities seen in timing 
data is free precession. This motion induces modulations in 
the spin down torque. Part of this modulation is a smooth 
function of precessional phase. However, part of the spin- 
down torque is either on or off, the on state corresponding 
to the energy threshold of some avalanche-like process being 
exceeded at some point on the stellar surface or magneto- 
sphere. 



4 COMPARISON WITH OBSERVATIONS 

Having constructed a simple free precession hypothesis, we 
now attempt to evaluate it in the light of recent pulsar ob- 
servations. We will draw our data set fro m the pulsar s whos e 
timing features w ere studied in detail in lLvne et al.l (J2010T ). 
ILvne et all (I2010T) report ed on 18 pulsars in total; 17 taken 
from iHobbs et all 1)20101 ) , and al so the intermittent p ulsar 
PSR B1934+21, first reported in iKramer et all l|2006l) . For 
each pulsar we computed the quantity 



N 



T 



(20) 



where T s is the time-span of the observations and P mo d the 
reported modulation period. The quantity N is therfore the 
number of modulation cycles that fit into the total data 
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span, and is (one) measure of how convincing the evidence 
is for harmonic structure in the timing data. Given that the 
data have been fit to timing models that include the pul- 
sar spin frequency v, frequency derivative z>, and (in some 
cases) second derivative i>, but no higher derivatives, some 
low order polynomial terms will inevitably remain in noisy 
data, even in the absence of precession, making the inter- 
pretation of low- TV oscillations as real features dangerous. 
We wi ll therefore make use of those pulsars from lLvne et al.l 
(J2010f ) that have N > 3; the exact cut-off of 3 is somewhat 
arbitrary. There are 15 such pulsars; we are not including 
pulsars B2035+36 (N = 0.37), J2043+2740 (N = 1.1) or 
B0950+08 (JV = 2.6) in our sample. We gi ve, in Table ED 



below , data for these 15 p ulsar s, taken fro m Krame r et al 
l|2006h . lHobbset al.1 (|2010h and lLvne etafl (|2010l 1. 

We will use these pulsars to examine our precessional 
hypothesis. 



4.1 Explaining the modulation period 

In the precession interpretation, we identify P mo( j with Pf p , 
the free precession period of a biaxial star, as described in 
section [31 modulo allowance for the fact that real neutrons 
star are not simple rigid bodies. When allowance is made for 
the fact that a real star consists of an elastic crust contain- 
ing a fluid core, the ratio of spin to free precession periods 
becomes 

P _ h-h 
Pi P Ic 



(21) 



where I3 — 7i refers to that part of the asymmetry in the 
moment of inertia tensor, generated by strains, that follows 
the precessional motion, while 7c refers to the crust plus all 
other parts of the star that are coupled to it on timescales 
shorter than the spin period (see I Jones fc Anderssonl (| 200 il l 
and references therein). 

The asymmetry 73 — 1\ we will continue to parameterise 
as a dimensionless number 

eai = ^ (22) 

where I* denotes the total stellar moment of inertia. The 
size of £31 will depend upon the strain mechanism at work, 
i.e. whether it is elastic or magnetic, as will be discussed 
below. 

The size of 7c requires some thought, as its value de- 
pends on the coupling between the crust and the rest of the 
star. If the entire star is threaded by a magnetic field, then 
the interior charged plasma (electrons and protons) will be 
tightly coupled to the crust ( Eassonl 1 19791 1. with a coupling 
timescale much shorter than any of the timescales of interest 
here. The crust and interior plasma can then be thought of 
as a single component, so that Ic at the very least consists of 
the moment of inertia of the crust, plus a small extra piece 
contributed by the interior plasma. 

It is then necessary to consider the interaction between 
the crust/plasma and the core neutron fluid. In the case 
that the core neutrons are normal (i.e. not superfluid) then 
they will be strong ly coupled to the charged plasma in a 
frictional-like way (|Bavm et al.lll969l ). The case of a two- 
component star, with its components coupled by a frictional 
torque proportional to t he vectorial angular v elocity differ- 
ence, was considered bv lBondi fe Gold! (1 19551 ). who showed 



than if this frictional coupling acts on a timescale of order 
the spin period P or less, the whole star, curst plus fluid 
core, will precess as one. Conversely, if the frictional cou- 
pling acts on timescales much greater than P, then the fluid 
core does not participate in the precession. It follows that in 
the case that the core neutrons are normal, the whole star, 
core plus crust, precesses as one, giving Ic ~ I*. 

On the other hand, if the interior neutrons are super- 
fluid, the interaction between crust and core is mediated by 
a process known as mutual friction. This can be provided 
by the ch arged plasm a scattering off the neutron rotational 
vortices (JAlpar fc Sa uls 1988). The precession of a two- 
component star, with the components cou pled by a mutual 
frictio n torque, was considered in detail bv lSedrakian et al.l 
(1999). In the limit of sufficiently weak mutual friction, a 
long period precessional mode still exists, with the neutron 
fluid decoupled from the precession, so that 7c would then 
be of order the crustal moment of inertia, i.e. Ic/I* ~ 10~ 2 
l|Sedrakian et al.ll 19991 ; iGlampedakis et alj|200il2009l ). 

However, it has been argued that there will be an even 
stronger interaction between the crust and fluid core if the 
interior protons are superconducting: then the magnetic flux 
tubes that make up the superconducting field will interact 
strongly with the neutron vortices. Link has argued that this 
interaction is so strong that there will be no long-lived long 
period free precession, forcing one to conclude that either 
the pulsar observations do not correspond to precession, or 
else the ne utron vorti c es do not coexist with the magnetic 
flux tubes ([Lin k 2003, 20061). In the latter case, one might 
conclude that the interior neutrons are normal, returning us 
to the Ic/I* ~ 1 case described above. Alternatively, one 
might conclude that the neutron vortices are spatially sepa- 
rate from the precessing crust and magnetic field, decoupling 
then from the precession, implying Ic/I * ~ 10~ 2 . Note that 
the recent cooling observations of Cas A (jHeinke fc Holl2010l ) 
suggest that both neutron superfluidity and proton super- 
condu ctivity are present in all but the youn gest of neutron 
stars (|Page et al.ll201ll ; IShternin et al.ll201ll ). favouring the 
second of these suggestions. 

Given the uncertainty in the correct value of Ic/I*, we 
will explicitly introduce the ratio Ic/I* in our calculations: 



P _ h-h 
Pi P 7 C 



hi* 



7 C =e31 lr 



7c 



(23) 



and consider both the cases Ic/I* ~ W~ and Ic/I* ~ 1 in 
our discussions. 

It remains to consider the value of the ellipticity param- 
eter e3i. For an unmagnetised fluid star £31 is zero; non-zero 
£31 values are generated only by strains, which may be ei- 
ther elastic or magnetic in origin. We will now consider the 
different possibilities. 

Let's first consider elastic strains, which will be sup- 
ported by the solid crust. For a star whose crust is relaxed 
when it has a shape parameterised by £0, the 'zero strain' 
or 'reference' oblateness, we expect 



£31 



be , 



(24) 



where fc is a measure of the rigidity of the crust 
llPines fc Shahaml [l972l : IJones fc Anderssonl l200ll ). In 
ICutler et al.l (|2003h it is estimated that b ~ 10" 7 for typ- 
ical neutron stars. 

For a star with a completely relaxed crust, its refer- 
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Table 1. Pulsar data used in our analysis. T he pulsar spin period P and modulation 
period P m o d are taken from l Lvne et al,l l|201Cf ). while the observation timespans T a are 
taken from Hobbs ct al. (2010), with the exception of PSR B1931+21 whose timespan was 
taken from Kramer ct al. (2006). The quantity N, defined in equation H20I I. is a measure 
of the significance of the reported modulations. The quantity T B AF is a measure of the 
stability of the modulation period, as discussed in secti on [4.21 AF is the uncertainty in 
the modulation frequency, taken from lLvne et al.l 1120101 ). 



ulsar name 


P (seconds) 


^mod (years) 


P/Pmod 


T s (years) 


N = n/P^ 


T S AF 


0740-28 


0.167 


0.370 


1.425e - 08 


20.60 


55.62 


4.12 


0919+06 


0.431 


1.613 


8.450e - 09 


27.00 


16.74 


1.08 


1540-06 


0.709 


4.167 


5.386e - 09 


19.70 


4.73 


0.39 


1642-03 


0.388 


3.846 


3.190e - 09 


35.70 


9.28 


2.50 


1714-34 


0.656 


3.846 


5.399e - 09 


15.20 


3.95 


0.61 


1818-04 


0.598 


9.091 


2.082e - 09 


35.10 


3.86 


0.35 


1822-09 


0.769 


2.500 


9.737e - 09 


19.70 


7.88 


1.38 


1826-17 


0.307 


3.030 


3.207e - 09 


18.20 


6.01 


0.36 


1828-11 


0.405 


1.370 


9.357e - 09 


18.70 


13.65 


0.37 


1839+09 


0.381 


1.000 


1.207e - 08 


25.30 


25.30 


3.79 


1903+07 


0.648 


2.778 


7.383e - 09 


19.60 


7.06 


2.55 


1907+00 


1.017 


6.667 


4.829e - 09 


23.60 


3.54 


0.47 


1929+20 


0.268 


1.695 


5.008e - 09 


17.20 


10.15 


0.34 


1931+24 


0.814 


0.076 


3.373e - 07 


1.67 


21.88 


1.17 


2148+63 


0.380 


3.030 


3.969e - 09 


19.80 


6.53 


1.39 



ence shape will be the shape of the equivalent rotating fluid, 
which has an oblateness of order the kinetic energy divided 
by the gravitational binding energy: 



6n 



MR 2 Q 2 
GM 2 /R 



2.1 x 10" 



(25) 



It follows that for a star whose crust is perfectly relaxed at 
its current rotation rate we can set eo = £n to give 



e 31 « ben « 2.1 x 10" 



10- 



which translates into a modulation period 

-T\ / p\ 3 

Is 



mod 



lyear 



1.5 x 10" 



10" 



(26) 



(27) 



In reality we would expect a spinning-down star to have 
a larger value of £31 than this, as the crust will be strained, 
retaining a memory of an earlier faster rotation rate. The 
exact picture depends rather subtly on the spin period at 
the moment when the crust first solidifies, and the breaking 
strain. To understand why, it is useful to consider the fol- 
lowing evolutionary picture. Suppose the crust first solidifies 
when the star is spinning fast. The crust is in fact believed to 
form very shortly after birth, so we will label the ellipticity 
of this relaxed star as eo = en, birth, with the understand- 
ing that, technically, 'birth' refers to the formation of the 
crust, not the slightly earlier birth of the star. As the star 
then spins down, a strain builds up; if the strain reaches the 
breaking value ttbrcak, a starquake occurs, resetting eo to a 
new value appropriate to the current rotation rate. This pro- 
cess can occur many times, such that eo ~ en throughout. 
However, when the star has spun down to the point that 
eo ~ Ubrcak, further spin-down will not generate sufficient 
strain to cause further starquakes, so that the strain is then 
'frozen in' at the level of Mbreak- Algebraically we have 



C3i ~ ben for P < P crit , 

631 ~ &«break for P > P cr it , 



(28) 
(29) 



where the critical dividing spin period P cr it is given by 
en (Pcrit) = ^brcak so that 



1.5 x 10" 



10- 1 



^brcak 



1/2 



(30) 



where our parameterisation of Ubrcak is motivated by 
the very high br e aking strains computed recently by 
iHorowitz fc Kadaul (|2009l ). 

However, if the breaking strain really is as high as this, 
the critical period P cr it of equation (|30[) is probably shorter 
than the actual birth spin period of pulsars, invalidating the 
above evolutionary scenario: spin-down would not induce 
sufficient strain to break the crust and so starquakes would 
not occur as the star spins down. Instead, the zero strain 
oblateness remains fixed at its birth value, eo = en, birth, for 
the entire stellar lifetime. In this case 



£31 



btn, birth 



(31) 



Once such a star has spun-down by a factor of a few, the 
level of strain would reach u « en, birth, below the breaking 

Strain Mbreak- 

Taking all this into account, we can parameterise the 
ellipticity: 



631 



631 



2.1 x 10" 



1.0 x 10" 



10- 7 
b 



for P< Pcrit, 
for P> Pcrit- 



10- 7 / vio- 1 

with corresponding modulation periods given by 



(32) 
(33) 



1 year 

tmod 

lyear 



1.5 x 10° 



10- 7 



3.2 



10" 



ai 

u 



for P < Pcrit, 


(34) 


for P > Pcrit- 


(35) 



Pulsar state switching, timing noise and free precession 7 



with the strain u in the P > P cr it case equal to 

U ~ Ubroak for Pbirth < -Pcrit, 
U ~ £n, birth for Pbirth > Pcrit- 



(36) 
(37) 



It follows that in the case of rapidly spinning star whose 
crusts can be approximated as relaxed we have P mo d oc P . 
For slow spinning stars with a fixed level of strain, P mo< j oc P, 
to the extent that all stars have identical breaking strains 
or identical birth spin periods. 

In the case of magnetic field deformation the situation 
is simpler. In the case of a magnetic field of typical strength 
B threading non-superconducting stellar fluid, the expected 
cllipticity is of order the magnetic binding energy to the 
gravitational binding energy: 



£31 

which leads to 

-*mod 



B 2 R A 
GM 2 /R 



1.9 x 10" 



B 



10 12 G 



1 year 



1.6 x 10 4 



10" 



(38) 



(39) 



\l/2 



where we have used B = 3.2 x 10 (PP ) '" (gaussian cgs 
units; see iLvne &: Graham-Smithl ((2006J)) and have param- 
eterised using values of B and spin-down rate P typical for 
the set of pulsars under consideration. 

In the case of superconductivity these results are modi- 
fied by the replacement of one factor of B with H c , the 'crit- 
ical field', believed to approximately equal to 1 15 G, inde- 
pend ent of the star's individual field strength B (|Bavm et al.l 
1 19691 ). The above results are then amended to 



£31 



BH C R A 

GM 2 /R 



1.9 x 10" 



B 



io i2 g; vio is g 



(40) 



This leads to 



16 



1/2 



10" 



1\ 1/2 



10 1S G 



1 year ^ 1 s J \ P J \ H c 

or, in terms of spin-down age, r sc i = P/(2P) 



lyear 



4.1 x 



T s d 



10 6 years 



1/2 



10 15 G 



(41) 



(42) 



(As mentioned previously, given that we are attempting to 
model long period free precession, the existence of super- 
conductivity is problematic, as a strong interaction between 
the magnetic flux tubes (not included in our model) that 
are believed to make up the magnetic field in a supercon- 
ducting star with a coexisting neutron superfluid vortex 
array would tend to increase the precession frequency to 
very high values, or else damp the precess ion very rapidly 
l| Jones fc Andersson|[200l1 : lLinkll2003l . l2006l ). A possible res- 
olution would be for the neutrons to be normal.) 

We therefore see that within the precessional hypothesis 
there are four different P mo d(P, P) relationships, depend- 
ing upon the mechanism at work: precession of an elasti- 
cally relaxed star (equation (|34[0 . precession of an elastically 
strained star (equation (1351) ). precession of a magnetised but 
non-superconducting star (equation (|39|) ) or precession of 
a magnetised superconducting star (equation (|42|) ~). These 
are not all mutually exclusive as both elastic and magnetic 
strains can contribute to the overall strain, while the two 
quoted elastic results represents the expected extremes of 



Figure 1. Plot of modulation period P mo d in years verses spin 
period P in seconds for 15 pulsars from Table 1 of Lvne et al.l 
J201Ch . Error bars taken from Table 1 of lLvne et ail J20ld ). The 
solid line is P/P mo d = 5 X 10 — 9 , and is not a fit to the data; 
a linear relation of this form is expected for the precession of 
stars with identical elastic strains generated by spin-down from 
an earlier faster rotation rate (see equation 1351 1. 




crustal strain, if one sets the strain to its breaking value 

^brcak- 

Clearly, the relaxed crust and normal magnetic field sce- 
narios can be rejected immediately, as the predicted modu- 
lations periods of equations (|34[) and (|39p are far in excess 
of these observed, regardless of where the ratio Ic/I* lies 
in the interval 10 -2 — > 1. In contrast, in the cases of free 
precession of a star with a strained crust or a superconduct- 
ing star, the periodicities of equations (|35p and (|42[) are of 
the right order of magnitude to explain the data, although 
in the latter case we would need Ic/I* ~ 1. To gain further 
insight we can make use of the differ ent scaling of P mo d with 
P and P, and exploit the fact that ILvne et al.l (|201dh have 
collected data on multiple pulsars. 

In Fig[T]we plot the modulation period P mo d verses spin 
period P for the pulsars given in Table Q] If the free preces- 
sion of a strained star is at work, the data points should 
form a line of unit gradient; to guide the eye we have plot- 
ted a line corresponding to P/P mo d = 5 x 10~ 9 (see equatio n 
(J35j) ) . Note the errors bars, taken from ILvne et al.l (|2010r) . 
measure the width of the power spectra peaks in the spin- 
down rate. All but two pulsars lie within a factor of 2 of the 
line P/P m od = 5 x 10~ 9 . The exceptions are pulsar B0740-28 
which lies a factor of 3 below the line, and the intermittent 
pulsar B1931+24, which lies a much larger factor of about 
70 below the line. 

In Fig [2] we plot the modulation period P mo d verses 
the pulsar age r s d = P/(2P). If free precession of a star 
deformed by superconducting magnetic field is at work, the 
data should form a line of gradient 1/2; to guide the eye, we 
have plotted such a line, corresponding to H c = 2 x 10 G 
(see equation (|42|l ) . There is slightly more scatter about the 
indicated line in this case, with four pulsars lying a factor of 
3 or more off the line. PSR B1931+24 is again a clear outlier, 
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Figure 2. Plot of modulation period P mo d i n years v erses pulsar 



ager = P/(2P) in years for 15 pulsars fr om Table 1 o f Lvne et al,l 
j201Ch . Error bars taken from Table 1 of lLvne et ail ||2010|) . The 
solid line is the expected scaling of P mot j for presession of super- 
conducting stars deformed by magentic strains, and is not a fit to 
the data. It was obtained by setting H c = 2 X 10 15 G in equation 
@2}. 
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lying a factor of 30 below the line. Also, the H c — 2 x 10 15 
G line appears to be somewhat steeper than the data. 

To make a more quantitative statement as to the ac- 
tual scaling of Pmod with P and P we have carried out a 
least squares fit to the data, excluding PSR B1931+24 due 
to its clear outlier status. The scatter of the points about 
a straight line is likely to be dominated not be the mea- 
surement errors in P mo d but by intrinsic differences between 
the stars themselves, e.g. differences in mass (which affects 
crust thickness and hence 6), differences in fracture history, 
or differences in internal magnetic field arrangement. We 
can see no easy way of estimating the sorts of scatter these 
effects are likely to produce in the modulation period. We 
will therefore not make use of any error estimates in cal- 
culating the best fit lines. Assuming a relation of the form 



aP P c , and fitting for a, b, c leads to 



(P/ls) 



(p/kt 



-0.14± 0.094 



lyear 



3.79 x 10" 9 (1±0.14) 



(43) 



This scaling is close to that of equation (|35[) . i.e. close to that 
expected for the free precession of a star deformed by elastic 
strains, with the level of strain approximately equal in all 
stars. Unfortunately our lack of knowledge of the intrinsic 
scatter in P mo d as a function of P and P prevents us from 
making any statement as to the statistical significance of 
equation (|43|l . All we can say is that of the four different 
deformation mechanisms proposed, the strained crust one is 
the most plausible. Note that if we had assumed from the 
outset that P mo d was only a function of P, and therefore 
carried out a fit of the form P mo d = aP b the result would 
have been rather similar: 



(44) 



4.8 x 10- 9 (1 ±0.13)' 



fit to that data, we can then use equation (|35[) to estimate 
the corresponding strain; obviously the actual braking strain 
must be at least as large as this: 

P/Pnod \ /10' 7 \ (I_C 

5 x 10- 9 M b \I* 



^break ^ U ~ 5-0 X 10 



(45) 



This lower bound on the breaking strain can then be trans- 
lated, using equation (|30f> . into an upper bound on the crit- 
ical spin period P cr it; the birth spin period must equal to or 
less than this: 

. ,, ,V2 / r „\ -1/2 /c „ 1n _ 9 \l/2 
Pbirth < Pcrit < 2.1mS 



10~ 7 



5 x 10" 



P/Pu 



(46) 

If the whole of the fluid interior participates in the pre- 
cession, so that Ic/I* ~ 1, we can therefore say Ubrcak > 
br;5.0x 10" 2 and P birth < 21 ms, or / birth > 490 Hz. This 
high level of strain is in fact consistent with the high break- 
ing st rains suggested by the work of iHorowitz fc Kadaul 
(2009). However, this birth spin period is rather short: sim- 
ple extrapolation of the Crab pulsar back to its historically 
known birth date suggests a much longer birth period. If 
instead only the crust participates in the precession, so that 



Ic/h 



10" 



we obtain 



^brcak 



> 



5.0 x 10~ 4 and 



Having seen that P/P mo d ~ 5 x 10 9 is a reasonable 



Pbirth < 0.021 s, or /birth > 49 Hz, a more plausible bound 
on the birth spin rate. 

In either case, the roughly uniform level of strain for all 
the precession candidates is interesting; the data points in 
Figure[T]are scattered by a factor ~ 2 about the best-fit line. 
One possible explanation for the uniformity is simply that 
this reflects the actual breaking strain of neutron star crusts. 
The larger value quoted above, u — 5 x 10 -2 , is close to the 
0.1 value that has been obtained in the molecular dynamics 
simulations of Horowitz fc Kadaul (|2009T ). As discussed in 
IChugunov fc Horowitz! ( 20101 ). when finite-temperature ef- 
fects are taken into account, the breaking strain becomes 
a decreasing function of time. Extrapulation to the long 
timescale relevant to precessing stars is dangerous, but a 
value of Ubrcak ~ 5 x 10 -4 seems sensible. The scatter by a 
factor of about 2 about the best fit line would then corre- 
spond to minor differences in crust thickness and 'geological' 
history between the different stars. Alternatively, the unifor- 
mity in strain may reflect similar spin periods at birth, with 
crust-only precession implying a birth spin of about 50 Hz 
and whole-body precession implying a higher birth spin fre- 
quency closer to 500 Hz. The scaling of en, birth with the 
square of the birth spin frequency would require a scatter 
of no more than ~ y2 about the average birth frequency, a 
tight clustering indeed. Given this, it seems more plausible 
that the observed strain values really are giving us informa- 
tion about a universal neutron star crustal breaking strain, 
but a definitive statement cannot be made. 

If we were to interpret the behaviour of PSR B1931+24 
as free precession, despite its outlier status, we obtain the 
following. For Ic/I* ~ 1 we obtain the non-sensical strain 
u — 3.3. For Ic/I* ~ 10~ we obtain u — 0.033, which cor- 
responds to a birth period of less than 2.6 ms, or a birth 
frequency greater than 390 Hz. So, in this case, sensible re- 
sults can only be obtained for crust-only precession. If PSR 
B1931+24 is indeed undergoing free precession, we can pro- 
vide no explanation as to why is has such an anomalously 
short precession period. 

To sum up, the relation between spin and modulation 
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periods is of roughly the form to be expected for free preces- 
sion in the regime where the crust retains a significant level 
of strain, built up from spin-down from a faster rotation 
rate, but PSR B1931+24 seems to have too short a modula- 
tion period to be described by the same model parameters 
as the rest of the data. 



Figure 3. P—P diagram for the pulsar population, with the sub- 
set from Table[T]identified. For comparison, we have also indicated 
those pulsa r s that have been observed to glitch. Data taken from 
iLvne et alj J2010h and the ATNF database. 



4.2 Stability of the modulation period 

Free precession should be a smooth regular motion, with 
a period set by the stellar structure, as described above. 
A neutron star has a very high moment of inertia, so one 
expects the underlying clock to be extremely accurate. De- 
partures from perfect periodicity in the modulations, if large 
enough, would make a precessional interpretation problem- 
atic. There are, however, two complicating factors that need 
to be taken into account. Firstly, the data are of finite du- 
ration, so some uncertainty in the extracted frequency is 
inevitable. Secondly, the long-term periodic modulation is 
presumably not the only source of timing noise in the data; 
if sufficiently large, the other sources of timing noise can 
make extraction of the periodic modulation difficult, and 
possibly make it appear less regular than is actually the 
case. Indeed, one such source of noise will be the stochastic 
nature of the magentospheric contribution to the spin-down 
torque. It follows that only by averaging over sufficiently 
many modulation periods would one expect to recover the 
underlying high-precision clock. This last point is in close 
analogy with the problem of measuring pulsar spin periods, 
as most pulsars display large pulse-to-pulse shape profile 
variations, and only by averaging over many can an average 
profile be constructed, then allowing the extraction of the 
highly stable spin period. 

To investigate the stability of the modulation, we can 
examine the power spect r a of t he spin-down rate as given in 
Figure S2 of ILvne et al.l (|2010f ) . The spectra for most stars 
consist of multiple-peaks, and all peaks are of finite width 
in frequency space. This in itself is not necessarily a sign 
of wandering in the modulation frequency: observation of 
a periodic system over a finite time span will generically 
produce such power spectra. Even in the case where the 
process under consideration is perfectly periodic, the width 
of these peaks is ~ T s -1 , where T B is the observation span. 
This is to be com pared with the actu al peak widths, AF, as 
given in Table 1 of lLvne et al.l (|2010l l. A convenient quantity 
to consider is then the ratio of these two frequency widths, 
T S AF; a value of T S AF significantly in excess of unity would 
imply that the finite duration of the data span is not the 
main factor determining the width of the frequency peaks. 
We give this quantity in Table [l] 

As is clear from the table, the width of the vast majority 
of peaks is about that one would expect if the finite dura- 
tion of the observations is the dominant source of frequency 
uncertainty. The only pulsars with T S AF > 3 are B0740- 
28, and B1839+09, whose peak widths are about four times 
T^ 1 . On the basis of these power spectra it seems reasonable 
to say that in all or most cases there is no clear signature 
of wandering in the modulation period, and the data seem 
to be consistent with a highly stable underlying clock. That 
this is true for P SR B1931+24 was implicit in the analysis of 
iRea et al.l (J2008J), who attempted to explain the behaviour of 
this system by invoking a companion object, a planet or low 




"Normal" pulsars 
ditchers 
Harmonic timing n 



Period / s 



mass star. Note, however, that Shabanovahas reported some 
variation in pulse phase for PSRs B0919+06 and B1642-03 
(Shabanova 2010J, |2009j). For now we will simply note that 
these pulsars do not have high values of T S AF (1.1 and 2.5, 
respectively) so by this measure at least are stable. We will 
comment on this further in section 15.1.21 when we discuss 
slow glitches. 



4.3 Excitation of harmonic timing noise 

What makes the subset of pulsars given in Table Q] special, 
i.e. why is it these pulsars, and not others, display interest- 
ing harmonic structure in their spin-down rates? To address 
this issue we plot a P-P diagram in Figure |3l where we 
plot 'normal' pulsars, and also the pulsars with harmonic 
timing residuals. For purposes of comparison, we have also 
identified those pulsars that have been observed to glitch. 

Indeed, one possible explanation for why only a small 
subset of pulsars display harmonic timing residuals might be 
related to glitches: perhaps the precession was excited by a 
recent glitch. In PSR B1828-11, where the periodicity stands 
out cleanly, we can see that there has been no significant de- 
cay in the long-term modulation over the observation span. 
This implies that the stars in our subset might have been 
set into precession decades ago, quite possibly before the 
first pulsar observations were taken. However, this simple 
glitch model is clearly not correct: as it apparent from the 
P-P diagram, the vast majority of pulsars that have been 
observed to glitch do not display harmonic timing structure. 
Clearly, it is not the case that a glitch automatically leads 
to harmonic timing structure. 

However, an alternative explanation suggests itself, 
which makes use of the switching behaviour that is present 
in at least some of the precession candidates. As explained 
above, precession with wobble angle 9 will, in the case of 
a smoothly-varying torque function, produce torque varia- 
tions of fractional order 6. However, for those special pulsars 
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which are delicately balanced between the magneto-spheric 
states, an additional fractional variation 8T/T is present, 
where ST is the difference in spin-down torques between the 
two magnetospheric states. If this difference is sufficiently 
large, the free precession will be much more apparent in the 
timing data of these stars than in those that are not deli- 
cately balanced between magnetospheric states: the switch- 
ing serves to magnify the torque variations. This suggests 
the following explanation for the apparent rarity of observed 
precession: those 15 or so candidates are special by virtue of 
being finely balanced between different states. The observed 
switching is then not an argument against free precession, 
but in fact the mechanism that magnifies its effects to make 
it apparent in the radio timing data. 

Of course, this model requires some non-zero level of 
precession to provide the clock to regulate the switching 
mechanism. Unfortunately, we see no way of estimating how 
small a level suffices to make the switching mechanism work. 
Certainly, if the balance between the two magnetospheric 
states is sufficiently delicate, the required precessional am- 
plitude may be sufficiently small so the low level of unde- 
tected precession that presumably lies buried in regular pul- 
sar timing noise might suffice. That timing noise should in- 
deed contain some precessional component is natural: what- 
ever mechanism is responsible for producing timing noise, 
it is likely to excite some small degree of free precession 
i|Cordeslll993r h This would only fail to be the case if the 
mechanism produced only perfectly axisymmetric variations 
in moment of inertia or torque fluctuations which only had 
a component along the rotation axis. Alternatively, one can 
invoke a recent glitch to provide the necessary precession. 
The lack of visible precession in the known glitching pul- 
sars is then understood by none of them being sufficiently 
delicately balanced between different magnetospheric states. 

It would then follow that the amount of timing noise 
should be greater in the precessing subset than for the other 
pulsars (or at least in those precessing stars where sharp 
changes in spin-down occur), as this subset has the conven- 
tional non-precessing timing noise contributions present in 
all pulsars, plus the small smooth precession torque compo- 
nent present (but undetectable) in all pulsars, plus the sharp 
switching component seen only in a subset. To investigate 
this hypothesis, we plot in Figure|3]the rms timing noise 02, 
computed after fitting for v an d its first two time deriva- 
tives, as taken from Table 1 of IHobbs et al.l (|20101 ). versus 
characteristic age. As can be seen readily by eye, the sub- 
set do have somewhat larger amounts of timing noise than 
other pulsars of similar ages, but only by a factor of a few 
providing modest but by no means compelling evidence in 
favour of this model. 



5 DISCUSSION 

5.1 Other possible explanations 

We have argued that the free precession hypothesis is not an 
al ternative to t he ma gnetospheric switching model advanced 
in lLvne et al.l i|2010f) ; the two are compatible, with the pre- 
cession phase-dependent accelerating potentials causing the 
switching from one magnetospheric state to the other. How- 
ever, several other models have been advanced to explain 



Figure 4. Plot of pulsa r rms timing noise 02 verses pulsar age 
P/(2P). Data taken from lHobbs et a.1.1 120101) . where this quantity 
is labelled 02- The 15 pulsars of Table[T]with interesting harmonic 
structure in their timing properties are identified. 
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harmonic timing structure. Most were developed before the 
switching behaviour was reported. We conrnrent on them 
briefly here. 



5.1.1 Tkachenko waves 



As was first pointed out bv iRudermanl (|1970al ). there is, in 
addition to free precession, another form of stellar oscilla- 
tion that might give modulation periods of order a year: 
Tkachenko oscillations in a rotating superfluid. These are 
oscillations in the superfluid vortex array. The oscillations 
would result in a slight variation in the stellar moment of in- 
ertia. By angular momentum conservation, this would man- 
ifest itself in timing variations. The Tkachenko wave speed 
is given by c\ — kQ/Stv, where k is the quantum of vorticity 
given in terms of the neutron mass m n : k = h/2m n . For an 
oscillation of wavelength A this leads to a mode period of 
order 



Pr ~ 1.4 years ( — 



1/2 



A 



10 6 cm 



(47) 



More detailed multi-fluid cal culations bv lHaskelll i|201l[ ) and 
iNoronha fc Sedrakianl (|2008i ) confirm the existence of such 
long period oscillations (for some portion of the microphys- 
ical parameter space at least), leading these authors to sug- 
gest Tkachenko oscillations as an explanation for the ob- 
serve timing variations, alternative to free precession. Also 
of interest is the model advanced by Popov (2008), where 
Tkachenko oscillations couple to free precession in those 
stars where Pf p « Pp. 

The attractive feature of an explanation based on 
Tkachenko oscillations lies in its lack of fine tuning — only 
basic properties of vortices go into equation (|47[) . However, 
there are a number of problems. Firstly, the best-fit scaling 
between P mo( j and P is much closer to P mo d oc P than to 
Pnod oc P 1 ' 2 , favouring precession over Tkachenk o waves; 
see equation (|43[1 . Secondly, the long periodicities of lHaskelll 
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J2011I ') and iNoronha fe Sedrakianl (|2008h are found by a 
plane wave analysis on a uniform background, and, crucially, 
for wave vectors exactly perpendicular to the vortex array. 
For wave vectors not perpendicular to the rotation, the os- 
cillation frequency is increased to a value of order the spin 
frequency. In a real star a global mode must be calculated in 
a spheroidal-like geometry, with regularity conditions at the 
centre and boundary conditions at the spheroidal surface. 
It is not clear if enforcing such a geometry will still allow 
for long-period oscillations. Finally, if the Tkachenko oscil- 
lations are to cause magnetospheric switching, presumably 
significant surface motions of the crust must occur, to create 
the time- varying surface electric fields that we have argued 
are needed. Such crustal motions will provoke elastic restor- 
ing forces, again raising the mode frequency to too high a 
level. 



ing data for the pulsars in question is not fatal for the pre- 
cession interpretation of their behaviour. 

Finally, what of the lack of pulse profile variation in PSR 
B0919+06? As first pointed out bv lRudermanl (ll970bT ). pre- 
cession with wobble angle 9 should result in the observer's 
cut through the beam varying by an angle ±9 once per pre- 
cession cycle. Why is this not apparent? A possible expla- 
nation might lie in magnetospheric switching. If switching 
occurs in this pulsar then, as explained above, the necessary 
wobble angle 9 might then be very small. If the component 
of the magnetospheric emission that is being switched on 
and off does not produce emission toward Earth, then only 
the smooth beam variation caused by the small precession 
angle would be apparent, which might be too small to be 
discernible. 



5.1.2 Slow glitches 

Shabanova has attributed the interesting timing behaviour 
of three of the pulsars in our sample (PSRS B0919+06, 
B1642-03 and B1822-09) to the phenomenon of 'slow 
glitches', where there is a slow rather than step-like increase 
in spin frequency, over a period of a few hundred days. As 
noted by Shabanova, there are a number of arguments in 
favour of the slow glitch interpretation over the precession 
one. Firstly, in PSR B1642-03, th e timing residuals are not 
quite periodic f Sha banoval 120091 ) . Instead the spacing be- 
tween peaks seems to vary by a factor of about 2.6 (but, in- 
terstingly, there seems to be some regularity in the this vari- 
ation; see below). Secondly, the residuals of PSR B0919+06 
do have a clear periodic structure, but seem to undergo a 
1/3 of a cycle shift during the observing period l|Shabanoval 
120101 '). Thirdly, in PSR B09 19+06 there is no observed varia- 
tion in p ulse profile, as wo uld be expected if precession were 
at work (|Shabanovall2010j ). 

While attractive, the slow glitch model has two draw- 
backs. Firstly there is no obvious mechanism to provide mag- 
netospheric switchi ng, something whi ch has been reported 
for PSR B1822-09 in lLvne et ail l^Old ). Secondly, there is no 
obvious mechanism to provide the well-defined periodicity 
(modulo the reported phase change) seen in PSR B0919+16 
or the approximately periodicity claimed for PSR B1642-03. 
Given this, we can ask just how serious are the difficulties 
with the precession interpretation given above. 

As regards the timing issue, as noted in section 14.21 a 
good measure of the timing stability is the product T S AF, 
as this averages over the duration span, minimising the ef- 
fects of the non-precessional contributions to the timing 
noise. The values of this quantity for the three pulsars are 
1.1, 2.5, 1.4, not significantly in excess of unity. So, averaged 
over the full data span, these objects are stable in modu- 
lation frequency. Also, for B1642-03, Shabnova notes that 
there is some regularity in the time intervals between peaks 
in the phase residuals, with a 60 year o verall envelope being 
tentatively identified (|Shabanoval l2009). That the residuals 
should have a more complex structure than a simple sinusoid 
is perfectly natural if the body has a significant degree of tri- 
axiality, as then the timing is sensitive to both the 73 — 7i 
and I2—I1 moment of in ertia differences l|Landau fc Lifshit a 
ll969l ; lAkgun et alJlJOOtlh . We therefore suggest that the tim- 



5.1.3 Asteroids and orbital companions 

Low mass binary companions are a natural and simple way 
of explaining long-term periodicities in pulsar data, but can 
a low-mass companion explain the magnetospheric switch- 
ing pheno menon? The b inary companion model was anal- 
ysed by IRea et al.l ([2008) for the intermittent pulsar PSR 
B1934+2. The companion was modelled as a low mass star 
on an eccentric ~ 35 day orbit. When at periastron the com- 
panion influences the pulsar emission, either by supplying 
additional plasma to power increased magnetospheric emis- 
sion, pushing the pulsar from the off to the on state, or by 
causing transitory wind accretion onto the pulsar, pushing 
it into th e propellor r egime , i.e. from the on to the off state. 
However, IRea et al.l (|2008r ) found that the orbit had to be 
extremely close to edge-on for this model to work, and that 
the likely wind from the low-mass companion was probably 
too low to affect the magnetosphere sufficiently. Whether 
or not the lower level of magnetospheric switching that has 
been seen in other pulsars can be explained by planetary 
companions is not clear. 

Regardless of this, the planetary hypothesis provides no 
good explanation for the distribution of modulation periods 
as a function of spin period of Figure [T] there is no reason to 
expect any particular relationship between the orbital and 
spin frequencies of companion and pul sar. This last point 
also a pplies to the asteroid model of ICordes fc Shannon! 
(2008), where ionisation of an infalling asteroid was used 
to supply the changing magnetospheric activity. 



5.1.4 Non-radial modes 

Recently iRosen et al.l (120111) have prop osed a model to ex- 
plain the findings of iLvne et al.l (J2010J ). Their model has a 
significant overlap with the one proposed here, as they too 
suppose that velocity variations in the surface of the star 
cause modulations in the magnetospheric spin-down torque. 
In their model precession is not the cause of this variation. 
Instead, the variations are due to the presence of non-radial 
oscillation modes with non-zero amplitudes at the stellar 
surface. The mode itself is of short period: the long-term 
variations are to be understood as being caused by varia- 
tions in the amplitude of excitation of a single mode, or the 
redistribution of energy over several modes. While attrac- 
tive, their model does not (as far as we are aware) provide 
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a natural explanation of why the mode excitation ampli- 
tudes would undergo the long-term fluctuations necessary 
to explain the long-term periodic structure seen in the tim- 
ing residuals or for the approximately linear scaling of this 
modulation period with the spin period (as per Figure [TJ. 

5.2 Further possible tests 

Having carried out a critical appraisal of the precession hy- 
pothesis, it is natural to ask what further tests might be 
carried out, that might strengthen or weaken the evidence 
in its favour. 

The most obvious test would be to identify more pre- 
cession candidates, particularly candidates wit h spin peri- 
ods ve ry different from those in the data set of iLvne et al.l 
(2010). The prospect of finding precess ion in long e r spin 
period pulsars does not seem good. The ILvne et al.l (120101 ) 
pulsars had periods in the range 0.1 < P < 1.0 second. As- 
suming that P m od ~ P/5 x 10 -9 , then a spin period P > 1 
seconds implies a modulation period P mo d > 6 years, too 
long to be visible in all but the longest of data sets. So, 
shorter spin period pulsars, with P < 0.1 seconds would 
provide more useful data. 

As was first pointed out by IShahaml (|1977f ). the inclu- 
sion of a pinned superfluid component alters the relation 
between P mo d and P drastically. If we assume small wobble 
angle 8 and that the moment of inertial difference 73 — 7i is 
much less than the moment in inertia of the pinned super- 
fluid, the P m od(P) relation becomes 



* m n d — -* 



Ic_ 

7sF 



*% 



7sF 

L 



(48) 



where 7sf and 7c denote the moments of inertia of the 
pinned superfluid and crust, re spectively (jSedrakian et al.l 
1 19991 ; iGlampedakis et all 120091 ). Note that by 'crust' we 
again mean the actual crust plus all other parts of the star 
coupled to it that participate in the precession. The ratio of 
Pmod to P then depends upon the extent of the pinning (i.e. 
on Tsf) and on how much of the star participates in the pre- 
cession (i.e. on 7c). There are several possibilities. In the case 
of pinning only in the inner crust, we have 7sf ~ 10~ 2 7*. 
If the core neutron fluid does not participate in the preces- 
sion, we have 7c ~ 10~ 2 , giving P mo d ~ P, while if the 
core neutron fluid does participate, we have 7c ~ I* giving 
Pmod ~ 10 2 P. However, if there is pinning in the core we 
have 7sf ~ 7* , so that if the core fluid does not participate 
(i.e. 7c ~ 10 _2 7*) we have P mo d ~ 10~ 2 P, while if the core 
fluid does participate (i.e. 7c ~ 7*) we have P mo d ~ P. 

Clearly, the modulation periods are very short in all of 
these cases as compared to the no superfluid pinning case, 
and may even be shorter than the spin period. It is not at all 
obvious how such a fast modulation would best be identified 
in pulsar data. What is clear is that observation of any one 
of these cases would then mean that an alternative non- 
precessional explanation must be found for the periodicities 
discussed in lLvne et al.l (|2010h . 

To summarise these ideas, and help guide observers as 
to which pulsar populations they might wish to target, we 
plot several curves of possible P m od(P) relations in Fig- 
ure [S] The line for a star precessing with an elastically 
relaxed crust corresponds to equation (|34fl with b — 10 -7 
and Ic/I* = 10 -2 . We also plot the precession period for a 



Figure 5. Some possible relations between modulation period 
P mo d and spin period P. The line labelled "Relaxed crust" was 
obtained from equation II27I I with b = 5 X 10 -7 . The line labelled 
"Stressed crust, no pinning" was obtained from equations 1341 1 
and gSJ with b = 1CT 7 , u = 5 X l(T 4 and Ic/I* = 10" 2 . The 
line labelled "Tkachenko" was obtained from 1471 with A = 10 6 
cm. The line labelled "Crustal pinning" was obtained by putting 
7sf/7c = 10 -2 in 1481 1. The line labelled "Core pinning" was 
obtained by putting Isf/Ic = 10 2 in 1481 1. 
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star which, at short rotation periods, can be approximated 
as having a relaxed crust, but at long rotation rates has 
a reference oblateness determined by the breaking strain 
Kbreak- These lines were obtained by setting b = 10 -7 , 
U — Ubrcak = 5 x 10~ 4 and Ic/I* = 10 -2 in equations (|34|l 
and (I35|l . The critical spin period dividing the two regimes 
is Pent ~ 20.5(5 x 10 _4 /u b rcak) 1/2 ms (see equation (|30j) ). 
As noted previously, for realistic breaking strains, the cur- 
rent precession candidates all lie on the P > P cr it part 
of the curve where P mo d is a linear function of P, with 
P/Pmod ~ kbreak (7* /7c) , so that observations of P mod 
and P only allow the extraction of the product of these 
factors. Data points corresponding to pulsars in the higher 
spin rate relaxed regime, where P mo d scales as P 3 , would 
allow the identification of the location of the knee, which 
would allow the value of Ubrcak to be extracted separately. 
Clearly, observations of more periodicities in more rapidly 
spinning pulsars would be needed to do this, and if Wbreak 
is as large as the recent microphysical modelling indicates 
l|Horowitz fc Kadaul 120091 1 . then few if any pulsars will re- 
side in this regime. Setting the moment of inertia ratio Ic/I* 
equal to 1 rather than 10 -2 would simply move these curves 
downwards to be a factor of 10 2 . The line for Tkachenko os- 
cillations is given by equation (I47|l . Its similar slope to the 
previously described line again suggests that observations of 
spin-down in milli-second pulsars (MSPs) would be very use- 
ful. However, the most striking feature of the plot is the well 
known many order-of-magnitude difference between the pin- 
ning and no-pinning predictions. Alternatively, for the pre- 
cession of stars deformed by superconducting fields, Figure 
[2] is the relevant one. In this case the shortest modulation 
periods would be ~ 0.1 years, for the youngest pulsars, while 
the longest modulation periods would be of order hundreds 



Pulsar state switching, timing noise and free precession 13 



of years, for the MSPs. Clearly, observers will need to be al- 
low for modulation over a wide range of timescales to cover 
all of the theoretical predictions. 

In addition to the use of pulse profiles and timing resid- 
uals, pulse polarisation is also a useful quantity to mea- 
sure when loo king for precession. According t o the rotating 
vector model l|Radhakrishnan fc CookeHl969l L because the 
pulse polarisation is tied to the magnetic field, which is in 
turn tied to the stellar crust, the precession would manifest 
itself as a longer-period variation superimposed on the po- 
larisation angle sweeps that occur once per rotation. Indeed, 
a r ecent study of p o larisa tion measurements was carried out 
bv lWeisberg et al.l (J2010|), who found that the position an- 
gle variations in 19 pulsars were significantly better fit by a 
sinusoid than by a constant. However, we have not carried 
out a statistical fit to this data as only one of the the 4 pul- 
sars in the 'Class I' i.e. the subset of pulsars considered most 
likely by the authors to display sinusoidal variations, meets 
our selection criterion of TV > 3. Clearly, analysis based on 
observations on a longer time baseline would be very inter- 
esting. 

It would be particularly interesting to see if the sharp 
change in pulse profile observed in some pulsars is also ac- 
companied by a sharp change in polarisation angle. It may 
be the case that the overall magnetic field structure is not 
significantly affected by whether or not the magnetospheric 
process is on or off, but rather simply acts as a guide, deter- 
mining the motion of the charged particles when the process 
is active. In this case, the polarisation angle may undergo 
perfectly smooth variations, tied to the underlying preces- 
sional motion, despite the sharp switches in pulse profile. 
However, it may be the case that the magnetospheric pro- 
cess, by switching on magnetospheric currents, does change 
the gross magnetic field structure. In this case sharp switches 
in polarisation would occur. Clearly there is some interest 
in looking at the polarisation variations in stars with sharp 
magnetospheric switches, but the presence of a sharp change 
in polarisation does not necessarily invalidate the preces- 
sional model proposed here. 



5.3 Comparison with gravitational wave 
spin-down upper limits 

In our simple biaxial precession model, an important quan- 
tity in determining the modulation period was the 73 — 7i 
asymmetry in the moment of inertia tensor; we defined 
631 = (73 — 7i)/7*, and this is related to the observable 
P/Pmod via equation (|23|l . so that 



P Ic 



(49) 



As was noted previously, a real star will be triaxial, not biax- 
ial, whic h would lead to a m ore complicated free precessional 
motion (|Akgun et al.ll2006T ). Nevertheless, the observed val- 
ues of the ratio P/P mo d are giving us information on the 
deformation of the stars' moment of inertia tensors, modulo 
the uncertainty in the value of the ratio 7c /I*. 

There is another aspect of pulsar physics that probes 
asymmetries in the moment on inertia tensor: gravitational 
wave emission. Specifically, the observed spin-down rates of 
non-precessing stars can be used to place upper limits on 



Figure 6. Histogram of estimated ellipticities 631 derived from 
precession candidates and, for comparison, histogram of e^x U P~ 
per limits derived assuming 100% gravitational wave spin-down 
for the known pulsar population. Two versions of the £31 his- 
togram are given: the filled one has Ic/I* = 10~ 2 , while the 
shaded one has Ic/I* = 1. The plot is truncated at ellipticities of 
10 -4 so as to restrict attention to the physically plausible values 
of eHf 1 corresponding to members of the MSP population. The 
non-MSP have unphysically large values of e!/ 1 "; the beginning of 
this distribution can be seen clearly in the rightmost bin. The 
outlier in the two £31 distributions is PSR B1931+24. 
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the 72 — 7i part of the inertia tensor. If we define 

h-h 

£21 = ; , 



(50) 



then an upper limit on this quantity is obtained by assuming 
that all of the kinetic energy loss associated with a spin- 
down rate P for a non-precessing pulsar with spin period P 
is going into gravitational wave energy emission, leading to 
the upper limit 



UL 
«21 



5c 5 PP J 



1/2 



32(2tt) 4 G7*J ' (51 ^ 

(see e.g. lAbbott et al.l (|2007l 0. This is an upper limit be- 
cause in reality some portion of the spin down, possibly 
most, is accounted for by electromagnetic rather than grav- 
itational wave emission. 

The ellipticity £21 can be non-zero only by virtue of 
strains within the star, so it is plausible that the mecha- 
nisms that limit the value of £31 (as discussed in section l4~T"|l 
will also limit the ellipticity £21. Given this, it is interesting 
to compare the estimates of £31 from precession, with the 
upper limits on £21 obtained from spin-down. A histogram 
of these quantities is given in Figure [S] We plot the his- 
togram only for ellipticities less that 10 -4 , for which the 
£2i L values correspond to members of the MSP population 
and are relatively plausible. (The non-MSP population part 
of the histogram consists of a second broad peak, extending 
all the way up to ellipticities of order unity, with a maximum 
at ellipticities of around 10 -2 . This is a consequence of the 
well known fact that the £2i L upper limit ellipticities for the 
non-MSP population are unphysically large.). 



14 Jones 



We plot two versions of the £31 histogram. In one we set 
Ic/I* = 10 -2 in equation (|49[) ; this peaks around values of 
a few times 10" 11 . In the other we set Ic/I* = 1 in equation 
(|49[) ; this peaks around values of a few times 10 -9 . Amus- 
ingly, this coincides with the peak in the histogram of upper 
limits on e^-t f° r the MSP population. On the basis of this 
seeming coincidence, one might speculate that Ic/I* really 
is close to unity and the physical mechanism limiting the two 
sorts of deformation is indeed one and the same. Given our 
analysis of section [4~T1 for the precession candidates we would 
then ascribe this to a highly strained neutron star crust, so 
that the product bubreak really is of order a few times 10 -9 , 
and that stars really do retain this level of non-axisymmetric 
deformation, so that the upper limits on £21 of the millisec- 
ond pulsar population are in fact estimates, not upper lim- 
its. If true, this would clearly be of interest for gravitational 
wave astronomers, and would add to the motivation of con- 
tinuing to target the MSPs in the advanced detector era. 
That neutron star crusts might be strong enough to sup- 
port such levels o f stres s was demonstrated prev iously by 
lUshomirskv et al.1 (|200Ch and ICutler et al.l (|2003T ); all that 
is new here is the additional evidence in support of this, 
in the form of the histogram of £31 values supplied by the 
multiple precession candidates. 

If it were to be the case that the spin-down upper limit 
cllipticities £2i L of the MSPs are indeed close to the actual el- 
lipticities, it would follow that electromagnetic energy losses 
are smaller then currently believed, which in turn would im- 
ply that the MSP magnetic field strengths are lower than 
the ~ 10 9 G inferred on the basis of magnetic dipole spin- 
down. Of course, the very fact that MSPs are visible as pul- 
sars at all implies that their magnetic fields are not zero, as 
some level of magnetic field is required to power the electro- 
magnetic pulsations. Furthermore, the X-ray and gamma- 
ray emission from MSPs is correlated with the total spin - 
down luminosity (|Possenti et alj I2002J : Icheng et aT1l2006h . 
This implies that whatever is causing the spin-down is also 
playing a role in producing the high energy electromagnetic 
emission, as would be expected for magnetic dipole spin- 
down, but not for gravitational wave spin-down from a star 
whose non-axisymmetry is supplied by elastic strain. This 
argues against a significant gravitational wave component 
to the spin-down. We simply note that the coincidence be- 
tween the peaks in the ellipticity histograms in the Ic/I* = 1 
case provides some small evidence that neutron star crusts 
might be strong enough to support the level of deformation 
needed to power gravitational wave emission from MSPs. 



6 SUMMARY AND CONCLUSION 

We have presented a model of pulsar timing where the har- 
monic structure of timing residuals seen in Lvnc ct al. (2010) 
is caused by free precession. We argued that the relationship 
between the modulation period and spin period was of the 
correct form to be explained by free precession of an elasti- 
cally strained star, the strain being due to the star retainin g 
a memory of a former faster spin rate (jCutler et al.l 12003). 
The details depend upon the fraction of the star that par- 
ticipates in the precession. In the case where only the crust 
participates (Ic/I* ~ 10 -2 ) we found the level of strain to 
be u ~ 5 x 10 -4 , and the stars must have been born spinning 



as least as rapidly as ~ 50 Hz. In the case where the whole 
star participates in the precession (Ic/I* ~ 1) we found 
the level of strain to be u at 5 x 10 -2 , and the stars must 
have been born spinning as least as rapidly as w 500 Hz. 
In either case, the roughly identical levels of strain in the 
different pulsars might then reflect the universal breaking 
strain that a crust can withstand. Alternatively, the nearly 
identical strain levels may reflect nearly identical spin rates 
at birth. 

We argued that the sharp switching in pulse profile that 
occurs in at least some of the candidate precessors is a con- 
sequence of the star being delicately balanced between two 
different magnetospheric states, with one state being statis- 
tically favoured over the other depending upon the preces- 
sion phase. That the link between magnetospheric state and 
precessional phase in only statistical was, we argued, con- 
sistent with the completely random mode-switching seen in 
some non-precessing stars. The sharpness of the switching 
reflects the energy of particles at some point in the magne- 
tosphere exceeding some critical threshold for an avalanche- 
like process, probably pair-production. 

Within this model, the magnetosphere variations, 
locked in phase with the precession, modulate the spin-down 
torque, serving to make free precession more readily de- 
tectable in such delicately balanced stars as compared to the 
rest of the population. It is therefore no surprise that magne- 
tospheric switching has been seen in precession candidates. 
The need for the star to be delicately balanced between two 
different magnetospheric states may play a role in explain- 
ing the rarity of the precession phenomenon in the pulsar 
population, although the precession excitation mechanism 
(perhaps glitching) may play a role too. 

We noted in passing that, in the case where the whole 
star participates in the precession, the required asymmetries 
in the moment of inertia tensor to explain the precession 
(£31 ~ 5 x 10~~ 9 ) are similar to the upper limits on asym- 
metry £21 that come from assuming the gravitational wave 
driven spin-down of the millisecond pulsar population. This 
may offer hope to gravitational wave astronomers that these 
upper limits may in fact be estimates, making the target- 
ing of the MSPs by future advanced detectors all the more 
interesting. However, if only the crust participated in the 
precession, the corresponding estimates of £31 are reduced 
by two orders of magnitude, comfortably below the upper 
limits on £21, removing this seeming coincidence. 

The analysis of crustal strain in precession candi- 
dates p resented here improves upon that of p revious works 
|jones fc Anderssonl l200ll ; ICutler et all 120031') by consider - 
ing a larger set of pulsars, as supplied bv lLvne et al.l (|2010T ) . 
and also by imposing the quality cut that at least three 
precession cycles should fit into the data set for the pre- 
cession claim to be taken seriously; we suggest that future 
investigators use a simil ar quality cut. Also , the analysis 
here improves on that of I Jones fc Anderssonl (|200ll ) by us- 
ing the_accurate value of the rigidity parameter b calculated 
bv lCutler et~ai] (|2003f) . 

There are a number of objections to the model presented 
here. Much of our argument is qualitative, rather than quan- 
titative. The only quantitative piece of evidence we have is 
the expected value of the modulation period for the slowly 
spinning stars discussed here, and its scaling with spin fre- 
quency, as illustrated in Figure [l] However, the spread in 
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spin periods (and modulation periods) was not large; more 
data points are needed to make this convincing. 

The re is also a serious th e oretic al prob lem. As pointed 
out b y I Jones fc Anderssonl (|200l[ ) and iLink fc Epstein! 
I|2001l ). the superfluid pinning that is believed to be responsi- 
ble for the larger Vela-like pulsar glitches would drastically 
reduce the modulation period, eliminating precession as a 
possible explanation of the o b served perio dic timing data. 
As has been argued by iLinkl (|2003l . 12000 ). if precession is 
indeed the mechanism responsible for the timing variations, 
it would be necessary for any superfluid parts of the star 
to be completely separated from the charged component, or 
at least that component that participates in the precession. 
This conflicts with microphysical modelling of neutron star 
interiors, which indicates that the interior superfluid is likely 
to coexist with a sold crust and, deeper within the star, with 
magnetic flux tubes. 

There is a caveat to be attached to this argument. The 
argument against the starquake model of Vela-like glitches 
is as follows. The strain relieved at each glitch is much larger 
than the strain build up by spin-down between glitches, 
so that the strain cannot possibly be replenished between 
glitches, i.e. a steady-state solution does not exist. However, 
the large crustal breaking strains suggested by recen t molec- 
ular dynamics solutions IjHorowitz fc Kadaul 12003 ) suggest 
a possible loop-hole. A typical Vela glitch releases a strain 



of order Av/v 



itr 
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once every 5 or so years, and is approximately 10 4 years 
old. This corresponds to a total of ~ 10 3 glitches over its 
lifetime, releasing a total strain of 10 -3 , two orders of mag- 
nitude smaller that the estimated breaking strain. So, the 
Vela and similar pulsars may simply be in the process of 
divesting themselves of crustal strain built up by spinning 
down from a much faster rotation rate that applied earlier 
in their lives. It follows that providing one is willing for the 
Vela and similar glitching pulsars to be in such a non-steady 
state situation, their large glitches can have a non-superfluid 
explanation. However, we can see no good reason why only 
a small fraction of the total strain should be relieved in each 
cracking event. 

To sum up, we feel that the precession model advanced 
here is plausible, and gives a sensible account of the ob- 
served timing data. However, the relatively modest number 
of pulsars with observations of periodic variation in spin- 
down, and their rather small spread in spin frequency, make 
forming a firm conclusion impossible. A larger sample of 
pulsars with harmonic timing residuals is needed, spanning 
a wider range of spin and modulation periods. By placing 
more points on the modulation period verses spin period 
plot (Figure 0), it should be possible to discriminate be- 
tween the precession-based model proposed here and the 
handful of others that have appeared in the literature. Also, 
a better picture of the stability of the clock that lies be- 
hind the timing data would be invaluable: modulo any irre- 
versible changes in stellar shape, precession should provide 
high frequency stability; a clear demonstration of a lack of 
modulation stability would argue against the precessional 
hypothesis. 
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